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ABSTRACT 

Measurements of the Cosmic Microwave Background Radiation (CMBR) provide a 
powerful tool for measuring the primary cosmological parameters. However, there is 
a large degree of parameter degeneracy in simultaneous measurements of the matter 
density, O^, and the Hubble parameter, Hq. In the present paper we use the presently 
available CMBR data together with measurements of the cosmological baryon-to- 
photon ratio, 77, from Big Bang nucleosynthesis, and the relative mass fraction of 
baryons in clusters to break the parameter degeneracy in measuring fi™ and Hq. We 
find that present data is inconsistent with the standard fJ = 1, matter dominated 
model. Our analysis favours a medium density universe with a rather low Hubble 
parameter. This is compatible with new measurements of type la supernovae, and the 
joint estimate of the two parameters is 57^ = 0.45lg'Q7 and Hq — 391^3 km s~^Mpc~^. 
We stress that the upper bound on the Hubble parameter is likely to be much more 
uncertain than indicated here, because of the limited number of free parameters in 
our analysis. 
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1 INTRODUCTION 



Barbosa 199^; ^ond and Jaffe 199^; [Bernardis et al.l997 



Tegmark 1998| ; [Lineweaver 1998| ; [White 1998[ )' 



Many recent measurements seem to indicate that our uni- 
verse is not a critical density, matter dominated Friedmann 
univers f Rather, the recent measurements nf type Ta s i ipia: 

novae _( Pcrl^^.ttor ^t nl 1 QQ«|; [Porl,T,„ttor ,>t nl 1 QQ' ' 



There is, however, a severe problem in that a change in 
one parameter can often be mimicked by suitable changes 
in a combination of other parameters ( Hu, Eisenstcin and 



tI hlarJ Tegmark 1999t [Eisenstein, Hu and Tegmark 1998a| ; [Eiscn- 



[navich U al 199^: [Riess e.t al 199^: [SrhmiHt et »/. 199^ stein, Hu and Tegmark 1998bj ). This is for instance very 



strongly suggest that, although the universe has a flat ge- 



ometry, the mat ter density is low (Perlmutter et al. 1998 



Riess et al. 1998). The energy de nsity is instead dominated 
by either a cosmological constant ( Carroll, Press and Turner 



1992), nr hy a gimrlar typp nf pnprgy ^untli nf^gathrn prpggiiT-p 



such as 



qiiintoss o nce fWang et al. 199^Zlatev, Wang and 



problematic when trying to simultaneously determine the 
cosmological matter density, Q.m,, and the Hubble parameter 
Hq. It has therefore been suggested that CMBR measure- 
ments should be combined with other constrai nts, coming 
for instance from large scale structure surveys ([Eisenstcin ' 



Steinhardt 1999|; [Caldwell, Dave and Stcinhardt 1998)) 



Hu and Tegmark 199 8a; Eisenstein, Hu and Tegmark 1998b; 
Webster et al. 1998|) or supernova type la rneasurements 



([Tegmark 199^ ; White 199S ; Lineweaver 199^; [Eisenstcin 



Nu pnerous investigations have shown that the fluctua- 



tion spectrum of the cosmic microwave background radia- 
tion (CMBR) provides a potentially very powerful tool for 
determining cosmological parameters (Jungman et al. 1996t 



HuandTcgmarkl998a ; Eisenstein, Hu and Tegmark 1998b; 



Bond, Efstathio u and Tegmark 1997 ; Zaklarriaga, Sperge] 
and Seljak 1997] ; |Hu, Eisenstcin and Tegmark 1999| ). Thus, 



one could hope that accurate measurements of these fluctu- 
ations could resolve the issue of whether or not the universe 



Efstathiou et al. 1999 ). In this paper we explore another 
possibility, namely how the cosmic baryon abundance can 
be used together with CMBR to provide tight constraints 
on fim and Hq. 

Galaxy clusters contain large amounts of hot. X-ray 
emitting gas from which the baryon mass to tot a l mass 
ratio can be measured (White and Fabian 1995; Pavid 



is dom nated by vacuum energy. Bas e d on present obser- 
yatio i is (Lineweaver and Barbosa 199^ ; pliveira-Costa et~al. 



Jones and Forman 1995|; Evrard, Metzler and Navarro 199C 



199£; Coble et al. 1999) there already exist a number of esti- 
mates of the cosmological parameters ([Efstathiou et al. 1999t 



Webster et al. 199E ; Hancock et al. 1998; Lineweaver and 



Evrard 1997 ). It turns out that measurements of this baryon 
cluster fraction can break the parameter degeneracy in- 
herent in the CMBR measurements, and although most 
of the present CMBR data is of relatively low accuracy 
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(JLineweaver and Barbosa 1998 Oliveira-Costa ei al. 



1999 



Coble ei al. 199£), it is still sufficiently accurate to pro- 



vide good constraints on the relevant cosmological models. 
We have used the presently available data, together with 
data from Big Bang nucleosynthesis and measurements of 
the cluster baryon fraction, and find that the observations 
are indeed strongly incompatible with a critical density mat- 
ter dominated model. Our results are easily compared with 
for inst ance the new supernova type la measurements ( Perl 



the CMBR likelihood contours are indeed narrowed signifi- 
cantly in the {Q,m, Hq) plane. The above way of breaking de- 
generacy will work when combining CMBR with any type of 
measurement that does not depend on Q,m and Hq as Q.mh^ ■ 
One such possibility is to use the cluster baryon frac- 
tion. It has long been known that measurements of the 
baryon cluster fraction favour a low density universe be- 
cause the measured fraction is so high that it ca nnot sup- 
port Hm = 1 without violating BBN co i istraints (Bludman 



mutter et al. 1998), and are found to be completely compat- 
ible. 

There still remains a quite large uncertainty on the de- 
termination of Q.m and Hq because of the low accuracy of 
the CMBR data. However, in the very near future, CMBR 
data of very high quality should become available from sev- 
eral different experiments. There is the balloon borne ex- 
periment BOOMERANG ( [Hanany 1997| ) which has already 
been fiown. Also, there are two new satellite experiments, 
MAP and PLANCK Q which wiU measure the fiuctuation 
spectru m very accurately on sub-degree scales. This should 



provide data which is accurate enough to diminish the un- 
certainties by an order of magnitude. 



2 BREAKING DEGENERACY 

As mentioned above, parameter extraction from the CMBR 
data suffers from some very large pa r ameter degeneracies 



( Eisonstein, Hu and Tegmark 1998a ; Eisenstein, Hu and 
Tegmark 1998bD . For instance it is not possible to con- 
strain Q.m a nd Hp separately, effectively only the combina- 
tion Q,mh^ (Eisenstein, Hu and Tegmark 1998a; Eisenstein 



Hu and Tegmark 19981) ; |Hu, Eisenstein and Tegmark 1999D 
where 



h = 



Ho 



100km s-iMpc" 



(1) 



can be measured accurately. 

However, as shown previously by several authors, it is 
possible to break this degeneracy by combining CMBR data 
with other, comp lementary, measurements. It was shown by 



Eisenstein et al. (Eisenstein, Hu and Tegmark 1998a| ; 



Eisen- 



atein, Hu and Tegmark 1998b[ ) that large scale structure sur 
veys like the Sloan Digital Sky Survey (SDSS) rl can probe 
the combination Om/i~^ . This means that the error el- 
lipses from such surveys are almost orthogonal to those of 
CMBR. The joint likelihood function should show a spike- 
like structure at the crossing of the two ellipses instead of 
the very elonga ted structure of either of the two individual 



measurements (Eisenstein, Hu and Tegmark 1998a; 



Eisen- 



atein, Hu and Tegmark 1998b| ). In fact, Webster et al. (Web 



3ter et al. 1998) have already performed a joint analysis of 
the present CMBR dat a together with da ta from the IRAS 
1.2 Jy galaxy survey (Fisher ei al. 1995). They find that 



* For information on these missions see the Internet pages for 

MAP ^ 

and 



( http: / /map.gsfc.nasa. gov 1 



(http://astro.cstcc.csa.nl/Planck/). 



PLANCK 



T See for instance [http: //www. astro. princcton.edu/BBOOK/ 



199^; IStcigman, Hata and Felten 1999| ; |White et^l. 1995| : 



Evrard 1997 ). A standard assumption in this game is to as- 
sume that the cluster baryon fraction is the same as the 
universal fraction, an assumption usually referred to as the 
fair sample hypothesis. Numerical simulation seem to justify 
this assumption. In fact recent simulations indicate that the 
cluster baryon fraction is slightly lower that that of the uni- 
verse as a whole (see (Steigman, Hata and Felten 1999) for 
a discussion). One problem is that observed clusters have 
diverging baryon fractions, d epending on their total mass. 
It is argued by Evrard ei al. (Evrard, Metzler and Navarre 



199C ) that this is nrost likely due to errors in the estimate of 
the total cluster mass, and that it can be corrected by use 
of statistical method s. In this pape r we shall use the results 
obtained by Evrard (Evrard 1997) for the universal cluster 
baryon fraction. 

The method for estimating f2m and Hq then works as 
follows: Big Bang nucleosynthesis (BBN) can be used t o 
measure the baryon-to-photon ratio (Kolb and Turner 1990), 
rj, which is related to fib by 

fibft^ = 3.66 X lO'^ 77. (2) 

Thus, nucleosynthesis only measures Qbh^ , not Qt- 

On the other h and, the cluster baryon fracti on, /s, is 
really measured as (Bludman 1998; Evrard 1997) 

r ,3/2 ^b ,3/2 /o\ 

fBh = —-h . (3) 

i 6777, 



By combining BBN measurements with the cluster baryon 
fraction we can therefore constrain the combination 



i irjr/t 



1/2 



(4) 



This combination of Qm and h is sufficiently different from 
Q.mh^ that when combining them it becomes possible to 
constrain both parameters well. 



3 MEASUREMENTS 

3.1 CMBR measurements 

In general, the fiuctuations in the CMBR is measured in 
ternrs of spherical harmonics 



T{e,4>)^^ai,^Yi^{e,4>), 



(5) 



where the coefficients are related to the power-spectrum Ci 
coefficients by 



Cl = {\air, 



(6) 



At present there is a host of different CMBR experiments, 
ranging from the largest scales (COBE), down to very small 
scales. The data that we us e are based on the compilatio n 
by Lineweaver and Barbosa ( Lineweaver and Barbosa 1998 ) , 
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but wit h the addition of the new Python V ( |Coble et al\ 
199S|) rpsults, and the resu lts from the QMAP experiment 
(pliveira-Costa et al. 199g|). 



3.2 BBN measurements 

The past few years have seen a very large fluctuation 
i n the estimate d baryon-to-photon ratio (see for instance 
(feteigman 199f ) for a review) . It has long been known that 



the primord ial value of deuterium p rovides a very sensitive 
probe of rj ( Kolb and Turner 1990 ), but the problem has 
been to measure the primordial deuterium abundance. In 
the local int erstellar medium , the abundance is quite well 



determined (Hata et al. 
D/H = 1.6 X 10"^ 



199! 



(7) 



but this can only really be used to provide a strict lower 
limit to the primordial abundance since deuterium is only 
destroyed, not produced, in astrophysical environments. 

Measurements of deuterium in quasar absorption sys- 
tems at high redshift have provided a completely new way of 
measuring the primordial abundance because such systems 
are chemically unevolved and should theref ore contain deu - 
terium abundances close to the primordial (Bteigman 1998). 

There have been two conflicting estimates of the deu- 
terium abundance in thes e systems, one which is much 
higher t h an the local value ( jongaila et al. 1994; [Carswell et 
' 1994|; iRugers and Hogan 1996| ; [Webb et al. 1997| ; [Tytler 
et al. 1998^ " 

D/H ~ 2 X 10"*, (8) 

and o ne which is only a factor o f two higher than the loc al 
value (iBurles and Tytler 1998a|; iBurles and Tytler 1998b|), 



D/H ~ 3 X 10" 



(9) 



There is growing evidence that the low deuterium value 
is the correct one, and as observational values we shall take 
the so-cal led Low-deuterium /High-helium data set, which is 
given by (JBurles and Tytler 1998a| ; [Buries and Tytler 1998b 
Izotov and Thuan 1998^ 



Yp = 0.245 ± 0.002 
D/H = (3.4 ± 0.3) X 10" 



(10) 
(11) 



This set of data is completely compatible with standard Big 
Bang nucleosynthesis for a baryon-to-photon ratio of 

-10 



r; = (5.1 ±0.3) x 10" 

or, in terms of baryon density Q,bh^ ~ 0.019 ± 0.001. 



3.3 Baryon cluster fraction 



(12) 



As mentioned above, Evrard (Evrard 1997) has calculated 
the universal cluster baryon fraction based on a large num- 
ber of clusters. His result is 



/s = (0.060 ± 0.003);i~^/^ 



(13) 



The la uncertainty is perhaps somewhat underes timated in 
the measurement, and we shall follow Hata et al. ( Steigman, 
Hata and Felten 1999) in assuming taking the la uncertainty 
to be twice as large, i.e. 0.006. The above value is derived 
from the gas fraction alone. The fraction in hot gas relative 



to coll apsed baryonic ob jects has been estimated by White 
et al. ([White et al. 1993[) to be 



M,. 



M, 



= 5.5/1-^/2 



gal 



(14) 



which is large enough to be insignificant. 



4 LIKELIHOOD ANALYSIS 

In order to estimate the underlying cosmological parame- 
ters, we have calculated a large number of synthetic mod- 
els which are then compared with the different data sets. 
Our synthetic models range through a large parameter space 
and have all been calculated using the publicly available 
CMBFAST code ( [Seljak and Zaldarriaga 1996| ). We re- 
strict the discussion to strictly flat models, i.e. models with 
Jim + ^A = 1. As free parameters we choose the normali- 
sation, Q, the matter density, flm, the baryon density, Qt, 
the Hubble parameter, h, and the spectral index, n. All the 
above parameters are allowed to vary with the restrictions 
described in table I. Altogether, 65000 independent CMBR 
spectra have been calculated, and, apart from the fact that 
we do not investigate open models, this analy sis is compa- 
rable to that of Lineweaver (Line weaver 199^). 

In a new analysis, Tegmark (Tegmark 1998) has calcu- 
lated the likelihood function in a larger, 9-dimensional pa- 
rameter space. The extra parameters are: the optical depth 
to reionisation, r, and the amplitude and spectral index of 
tensor fluctuations. Here, it was shown that relaxing the 
assumption about zero curvature significantly loosens the 
constraint on the combination Qmh^- In fact Tegmark finds 
that there is no upper bound on Q,mh? , even at the 68% 
confidence limit, contrary to our findings below. The reason 
for this is that the constraint relies to a certain extent on 
the amplitude differences in the power spectrum at different 
I. The information stored in this difference is sensitive to 
changing the spatial curvature, and therefore relaxing the 
flatness assumption leads to a much looser upper bound on 
Q,m.h? . This being said, there are good reasons for neglecting 
spatial curvature. Flatness is a generic prediction of almost 
all inflationary models, and furthermore the new data on 
Supernovae of type la c learly favour a spatially flat universe 
( [Perlmutter et al. 199^ ). 

From the same prejudice that has led us to consider only 
flat models, we have neglected possible tensor fluctuations. 
Almost all inflationary models predict only very small tensor 
fluctuations, so this is likely to be a good assumption. We 
have also neglected reionisation in our calculations. There 
is no justification for this since we know for a fact that the 
universe was reionised at fairly high redshift. As shown by 
Tegmark (Tegmark 1998), adding reionisation to the model 



parameters broadens the likelihood function, but does not 
move the maxima. Therefore we stress that our likelihood 
estimates should be considered optimistic and that they will 
broaden if reionisation is taken into account, but that they 
should still home in on the correct central values. 

In order to compare with the different data sets we then 
perform a x^ analysis. In general, x^ is given as 



JV 

/ 2\ i,ohs -^i, theory J 

i — 1 



(15) 
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Table 1. The free parameters used in our analysis 



Parameter 


range 


no. grid points 


Q2 


5-40 ^lK 


200 


n™ 


0.2-1 


17 


C6/l2 


0.002-0.030 


15 


h 


0.30-0.75 


16 


n 


0.7-1.3 


16 




70 
60 
50 
40 
0.2 



CMBR 



BBN + Clusters i 




1.0 



Figure 1. The likelihood function plotted in the {Q,m,Ho) plane 
after maximising over all other free parameters. The contours 
shown are 68%, 95% and 99% confidence limits. 



fidence limits correspond to Ax^ ~ 2.29, 6.18 and 9.21 re- 
spectively. The top panel shows the likelihood contours for 
CMBR data alone. As expected it shows a very elongated 
structure, roughly on the Q,m (x h~^ line. The middle panel 
shows how the likelihood contours are affected if one com- 
bines CMBR with the BBN estimate of Q.th'^. Clearly, there 
is only a quite small effect on the determination of flm and 
Ho, because Qb is not very degenerate with either of these 
two parameters. The real improvement only comes when one 
takes into account the cluster baryon fraction measurements. 
Doing this narrows the likelihood contours to a more vertical 
and much tighter structure. Intriguingly, the combined like- 
lihood function from this method favours a medium density 
universe with Sim around 0.5 and a very low value of the 
Hubble parameter. This is strongly incompatible with the 
standard Qrn ~ 1 CDM model, unless the Hubble param- 
eter is much lower than 0.3. However, as mentioned above 
our likelihood contours are somewhat optimistic, especially 
the upper bound on the Hubble parameter is sensitive to 
adding extra parameters to the analysis. 

It is also of interest to compare our results with the 
recent results from measurements of supernova type la's. 
The best available data come from the supernova cosmol- 
ogy project, who have m easured 42 high redshift supernovae 
( Perlmutter et al. 1998| ). From these they obtain stringent 
limits on Qm in a flat universe, but no constraint on the 
Hubble parameter. Their best fi t value for n^ in a flat uni- 
verse is (Perlmutter et al. 1998) 



i Lj] 



0.2810^08 (statistical) _o^o4 (systematic) 



(17) 



These data are compatible with our results, even at the la 
level, but only for relatively low values of the Hubble pa- 
rameter. Joint maximum likelihood for our analysis and the 
supernova data are shown in Fig. 2. The likelihood function 
has a maximum at 



(!:!™,/i) = (0.45,0.39), 



(18) 



but with a fairly large uncertainty on the Hubble parameter. 
At the 95% confidence level the Hubble parameter is only 
constrained to be less than 0.64 and, again, we stress that 
adding more parameters to the analysis could weaken this 
bound further. 



where F for CMBR data is equal to AT for different values 
of I. For BBN data, F is equal to either the helium or the 
deuterium abundance, ai is the uncertainty in the given data 
point. 

Here it has been assumed that the experimental errors 
are of Gaussian nature, which in reality they are not. There- 
fore the confidence limits that we obtain are not strict in a 
mathematical sense, but should rather be seen as indicative. 

If the measurements are furthermore assumed to be 
independent then the likelihood function is then given by 
L (X e ^ ' , and 



-^joint 



= Lc 



MBR^BBN.tJclustera. 



(16) 



In Fig. 1 we show the likelihood contours taking into ac- 
count different measurements. The likelihood contours have 
been based on Ax^, such that the 68%, 95%, and 99% con- 



5 CONCLUSIONS 

We have calculated the possible constraints on Q.m and Hq 
from combining CMBR data with data on the cosmic baryon 
abundance. It was found that these measurements are incon- 
sistent with the standard flat CDM model, even for very low 
values of the Hubble parameter. Our favoured region of pa- 
rameter space is consistent with re cent measurements from 
type la supernovae at the la level (Perlmutter et al. 1998). 
The joint likelihood function suggests that Q.m is close to 
0.4 and that a rather low value of the Hubble parameter is 
favoured. It is int eresting to compare our results with those 
of Webster et al. ( Webster et al. 1998| ), who did a j oint anal- 
ysis of t he CMBR data and the IRAS 1.2 Jy data ([Fisher ei 



al. 1991:) on large scale structure. They fou nd a best flt for 



the joint likelihood of ( Webster et al. 1995 ) 

(ilm, /l)cMBR+IRAS ~ (0.39 



^tri;^,o.53 



+0.05 \ 
O.14J 



(19) 
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70 


- 


\ CMBR i 






^ \ + BBN I 


60 


— 




\ \ + Clusters - 






\ \ + SN la I 


50 


- \ 


iN\ 1 


40 


f \k\\x i 



0.2 



0.4 



0.6 



.0 



Figure 2. Joint likelihood fun ction from combining the results 
of Rcf. (Pcrlmutter et al. 1998) with those of our analysis. The 
contours are 68%, 95% and 99% confidence limits. 



which is completely compatible with our estimate of 

V^im 5 /ijprcscnt analysis ^ V^-'^'^ — 0.07 ' ^•'^'^ — 0.13 / " \^^) 

In conclusion what we have shown in the present paper 
is that if the cosmic baryons are used in unison with CMBR 
data to constrain the cosmological parameters Q.m and h, 
the results are very close to other recent measurements us- 
ing different methods. This suggests that there is a con- 
vergence towards the true, underlying, cosmological model, 
which is apparently a flat, low density universe. The fact that 
so many completely independent methods all yield roughly 
consistent estimates make it increasingly unlikely that the 
data are contaminated by some large, unknown, systematic 
error. Of course when the new high precision data become 
available from the BOOMERANG, MAP and PLANCK ex- 
periments, these uncertainties should be resolved. 
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